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Abstract 

Assuming that the time-evolution of the self-consistent mean field is determined 
by five pairs of collective coordinate and collective momentum, we microscopically 
derive the collective Hamiltonian for low-frequency quadrupole modes of excita¬ 
tion. We show that the five-dimensional collective Schrodinger equation is capa¬ 
ble of describing large-amplitude quadrupole shape dynamics seen as shape coex¬ 
istence/mixing phenomena. We focus on basic ideas and recent advances of the 
approaches based on the time-dependent mean-field theory, but relations to other 
time-independent approaches are also briefly discussed. 


1 Introduction 

In this paper, we focus on low-frequency quadrupole motions which play the major 
role in low-energy excitation spectra. As is well known, various giant resonances 
appearing in highly excited states are well described by the random-phase approx¬ 
imation (RPA), which is a small-amplitude approximation of the time-dependent 
Hartree-Fock (TDHF) theory. In contrast to the giant resonances, low-frequency 
quadrupole vibrations exhibit characteristic features associated with superfluidity 
of the finite quantum system (nucleus), that is, pairing correlations and varying 
shell structure of the self-consistent mean field play essential roles mmm®- 

Quantum shape fluctuation, shape mixing and shape coexistence 

The low-frequency quadrupole vibrations can be regarded as soft modes of the 
quantum phase transition towards equilibrium deformations of the mean field. As is 
well known, in nuclei situated in the transitional region from spherical to deformed, 
amplitudes of quantum shape fluctuation about the equilibrium remarkably increase. 
This is the case also for weakly deformed nuclei where the gain in binding energies 
due to the symmetry breaking is comparable in magnitude to the vibrational zero- 
point energies. The transitional region is wide and those nuclei exhibit quite rich 


1 


excitation spectra. This is a characteristic feature of finite quantum systems and 
provides an invaluable opportunity to investigate the process of the quantum phase 
transition through analysis of quantum spectra. To describe such large-amplitude 
collective motions (LACM), we need to go beyond the small-amplitude approxima¬ 
tion (quasiparticle RPA) of the time-dependent Hartree-Fock-Bogoliubov (TDHFB) 
theory for superfluid systems. It is required to develop a microscopic theory of 
LACM capable of describing the varying shell structure associated with the time- 
dependent mean field with superfluidity. 

The spherical shell structure gradually changes with the growth of deformation 
and generates ‘deformed shell structures’ and ‘deformed magic numbers,’ that sta¬ 
bilize certain deformed shapes of the mean field. When a few local minima of the 
mean field with different shapes appear in the same energy region, LACM tunnel¬ 
ing through potential barriers and extending between local minima may take place. 
These phenomena may be regarded as a kind of macroscopic quantum tunneling. 
Note that the barriers are not generated by external fields but self-consistently 
generated as a consequence of quantum dynamics of the many-body system under 
consideration. Quantum spectra of low-energy excitation that needs such concepts 
have been observed in almost all regions of the nuclear chart [5]. When different 
kinds of quantum eigenstates associated with different shapes coexist in the same 
energy region, we may call them ‘shape coexistence phenomena .’ This is the case 
when shape mixing due to tunneling motion is weak and collective wave functions 
retain their localizations about different equilibrium shapes. On the other hand, if 
the shape mixing is strong, large-amplitude shape fluctuations (delocalization of the 
collective wave functions) extending to different local minima may occur. 

Collective rotations restoring broken symmetries 

As is well known, the central concept of the BCS theory of superconductivity 
is spontaneous breaking of the gauge symmetry and emergence of collective modes. 
The massless modes restoring the broken symmetry are called Anderson-Nambu- 
Goldstone (ANG) modes [6j El E] - As emphasized by Bohr and Mottelson, nuclear 
rotation can be regarded as an ANG mode restoring the broken rotational symmetry 
in real space HUH- 

In finite quantum systems such as nuclei, the rotational ANG modes may couple 
rather strongly with quantum shape-fluctuation modes. For instance, even when 
the self-consistent mean field acquires a deep local minimum at a finite value of 
/3, the nucleus may exhibit a large-amplitude shape fluctuation in the 7 degree of 
freedom, if the deformation potential is flat in this direction. Here, as usual, (3 
and 7 represent the magnitudes of axially symmetric and asymmetric quadrupole 
deformations, respectively. Such a situation is widely observed in experiments and 
called 7 -soft. Although the quantum-mechanical collective rotation is forbidden 
about the symmetry axis, the rotational degrees of freedoms about three principal 
axes are all activated once the axial symmetry is dynamically broken due to the 
quantum shape fluctuation. Rotational spectra in such 7 -soft nuclei do not exhibit 
a simple 7(7 + 1) pattern. Such an interplay of the shape-fluctuation and rotational 
modes may be regarded as a characteristic feature of finite quantum systems and 
provides an invaluable opportunity to investigate the process of the quantum phase 
transition through analysis of quantum spectra. 
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Thus, we need to treat the two kinds of collective variables, i.e ., those associated 
with the symmetry-restoring ANG modes and those for quantum shape fluctuations, 
in a unified manner to describe low-energy excitation spectra of nuclei. 

Five-dimensional quadrupole collective Hamiltonian 


Vibrational and rotational motions of the nucleus can be described as time- 
evolution of a self-consistent mean field. This is the basic idea underlying the unified 
model of Bohr and Mottelson mm ■ In this approach, the five-dimensional (5D) 
collective Hamiltonian describing the quadrupole vibrational and rotational motions 
is given by mm 


T v ib 


Hcoll — ^vib + ?rot + V ((3, 7 ), 

7)/3 2 + D Pl (p, 7)/?7 + ^T> 77 (/3, qjf 2 , 


Trot ~ 


I 2 


2J4(/?,7)' 


(1) 

( 2 ) 

(3) 


Here, (3 and 7 are treated as dynamical variables, and (3 and 7 represent their 
time-derivatives. They are related to expectation values of the quadruple operators 
(with respect to the time-dependent mean-field states) and their variations in time. 
The quantities (Dg^, D i g 7 , and D 77 ) appearing in the kinetic energies of vibrational 
motion, T v ib, represent inertial masses of the vibrational motion. They are functions 
of /3 and 7 . The quantities Ik and HkiPi 7 ) in the rotational energy T rot represent 
the three components of the angular momentum and the corresponding moments of 
inertia, respectively. Note that they are defined with respect to the principal axes 
of the body-fixed (intrinsic) frame that is attached to the instantaneous shape of 
the time-dependent mean-field. 

In the case that the potential energy V ((3, 7 ) has a deep minimum at a finite value 
of /3 and 7 = 0° (or 7 = 60°), a regular rotational spectrum with the 1(1+1) pattern 
may appear. In addition to the ground band, we can expect the (3 and 7 bands to 
appear, where vibrational quanta with respect to the f3 and 7 degrees of freedom 
are excited. Detailed investigations on the 7 -vibrational bands over many nuclei 
have revealed, however, that they usually exhibit significant anharmonicities (non- 
linearlities) [13]. Also for the /3-vibrational bands, it has been known mmm 
that they couple, sometimes very strongly, with pairing-vibrational modes (associ¬ 
ated with fluctuations of the pairing gap). Recent experimental data indicate the 
need for a radical review of their characters [5]. 


Collective quantization of time-dependent mean fields 


States vectors of time-dependent mean field are kinds of generalized coherent 
states, and we can rigorously formulate the TDHFB as a theory of classical Hamil¬ 
tonian dynamical system of large dimension [is moil [ 20 ]. Because time-evolution 
of the mean field is determined by the classical Hamilton equations, we cannot de¬ 
scribe, within the framework of the TDHFB, quantum spectra of low-lying states 
and macroscopic quantum tunneling phenomena such as spontaneous fissions and 
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subbarrier fusions. To describe these genuine quantum phenomena, we need to in¬ 
troduce a few collective variables determining the time-evolution of the mean field 
and quantize them. We refer this procedure ‘collective quantization.’’ 

For small-amplitude vibrations about an HFB equilibrium, it is well known that 
we can introduce collective variables in a microscopic way by solving the quasipar¬ 
ticle RPA (QRPA) equations. Single-particle spectra for a mean-field of a finite 
quantum system have rich shell structures, and thereby a variety of collective vibra¬ 
tional modes emerge. Even within the isoscalar quadrupole vibrations, two collective 
modes appear exhibiting quite different characteristics; the low- (usually first ex¬ 
cited 2 + ) and high-frequency (giant resonance) modes. One of the merits of the 
QRPA is that we can determine the microscopic structures of the collective coor¬ 
dinates and momenta starting from a huge number of microscopic (particle-hole, 
particle-particle, and hole-hole) degrees of freedom. We can thus learn how collec¬ 
tive vibrations are generated as coherent superpositions of many two-quasiparticle 
excitations. Examining the microscopic structure of the low-frequency quadrupole 
vibrations, we see that the weights of two-quasiparticle excitations near the Fermi 
surface are much larger than those in the mass quadrupole operators (see, e.g., 
Ref. (STj). This clearly indicates the importance of describing collective modes in 
a microscopic way. Another merit of the QRPA is that it yields the ANG modes 
and their collective masses (inertial functions). In this way, we can restore the 
symmetries broken by the mean-field approximation [8; 22) [23]. 

It has been one of the longstanding fundamental subjects in nuclear structure 
physics to construct a microscopic theory of LACM by extending the QRPA concepts 
to non-equilibrium states [Ml EH [ 26 ]. Below we shall briefly review various ideas 
proposed up to now for this aim. 


2 Basic ideas of large-amplitude collective mo¬ 
tion 

During the attempts to construct a microscopic theory of LACM since the latter 
half of the 1970s, significant progress has been achieved in the fundamental con¬ 
cepts of collective motion. Especially important is the recognition that microscopic 
derivation of the collective Hamiltonian is equivalent to extraction of a collective 
submanifold embedded in the TDHFB phase space, which is approximately decou¬ 
pled from other non-collective degrees of freedom. From this point of view we can 
say that collective variables are nothing but local canonical variables which can be 
flexibly chosen on this submanifold. Below we review recent developments achieved 
on the basis of such concepts. 

Extraction of collective submanifold 

Attempts to formulate a LACM theory without assuming adiabaticity of large- 
amplitude collective motion were initiated by Rowe and Bassermann m and Maru- 
rnori [28], and led to the formulation of the self-consistent collective coordinate 
(SCC) method [29]. In these approaches, collective coordinates and collective mo¬ 
menta are treated on the same footing. In the SCC method, basic equations deter¬ 
mining the collective submanifold are derived by requiring maximal decoupling of 
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the collective motion of interest and other non-collective degrees of freedom. The 
collective submanifold is a geometrical object that is invariant with respect to the 
choice of the coordinate system whereas the collective coordinates depend on it. 
The idea of coordinate-independent theory of collective motion was developed also 
by Rowe [ 30 j and Yamamura and Kuriyama [ 19 ] . This idea gave a deep impact on 
the fundamental question ‘what are the collective variables.' The SCC method was 
first formulated for TDHF, but later extended to TDHFB for describing nuclei with 
superfluidity Ell¬ 
in the SCC method, under the assumption that time evolution is governed 
by a few collective coordinates q = (qi, q-2 ,..., (?/) and collective momenta p = 
{PhP2, ■—,Pf ), the TDHFB states vectors are written as 

\m = \Hq,p)) = e ld{q ' p) \M, ( 4 ) 


or equivalently, 

where | <7^0) denotes the HFB ground state and 77 = (771,772, ■■■iVf) with 


Vi = Qi + iPi ), 


= -j=(Qi-iPi)- 


( 5 ) 

( 6 ) 


The TDHFB states |</>(f)) are required to fulfill the canonical variable conditions 

(HV,V*)\-^\^(V,V*)) = \v t, ( 7 ) 

I <!>{V,V*)) = ~\vu (8) 

which guarantee that ( q,p ) are canonical conjugate pairs. The one-body operator 
(7(77,77*) is determined by the time-dependent variational principle 

<5 ( 0 ( 77 , 77 *)| — H 1 ^( 77 , 77 *)) = 0 (9) 

with 

5 = + 00 ) 

and the canonical variable conditions. 

Making a power-series expansion of G with respect to (77,77*), 


G( 77,77*) = E E . (U) 

ij miTlj 

and requiring that the time-dependent variational principle holds for every power of 
(77,77*), we can successively determine the one-body operator This method 

of solution is called the (77,77*) expansion method. The collective Hamiltonian is 
defined by the expectation value of the microscopic Hamiltonian with respect to 
1^(77,77*)). Because (77,77*) are canonical variables, they are replaced by boson op¬ 
erators after canonical quantization. The lowest linear order corresponds to the 
QRPA. Accordingly, the collective variables (77*, 77*) correspond to specific QRPA 
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modes in the small-amplitude limit. It is important to note, however, that the 
microscopic structure of G changes as a function of (rp. rf *) due to the mode-mode 
coupling effects among different QRPA modes in the higher order. In this sense, 
the SCC method may be regarded as a dynamical extension of the boson expansion 
method [32] . Thus, the SCC is a powerful method of treating anharnronic effects 
to the QRPA vibrations originating from mode-mode couplings, as shown in its 
application to the two-phonon states of anharmonic 7 vibration [13L 33] . The SCC 
method was also used for derivation of the 5D quadrupole collective Hamiltonian 
and analysis of the quantum phase transition from spherical to deformed shapes 
[3T. and for constructing diabatic representation in the rotating shell model [35] ]. 

Solution with adiabatic expansion 

The ( 77 , rj*) expansion about an HFB equilibrium point is not suitable for treating 
situations such as shape coexistence, where a few local minima energetically compete 
in the HFB potential energy surface. For describing adiabatic LACM extending over 
different HFB local minima, a new method has been proposed [36]. In this method, 
the basic equations of the SCC method are solved by an expansion with respect 
to the collective momenta. It is called ‘adiabatic SCC (ASCC) method.’ Similar 
methods have been developed also by Klein et al. m, and Almehed and Walet 

[SSI- 

Let us assume that the TDHFB state vector can be written as 

\Hq,p))=e i ^ Pi ° i{q) m)). (12) 

Here, Q l {q) are one-body operators, called infinitesimal generators, and \4>{q)) is an 
intrinsic state at the collective coordinate q , called a moving frame HFB state. 

We determine the microscopic structures of Q l (q ) and \4>(q)) by the time-dependent 
variational principle 

H\(j){q,p)) =0. (13) 

Making power-series expansions with respect to the collective momenta p and 
retaining terms up to the second order in p, we obtain 

moving-frame HFB equation 


6 (cj)(q) | H M {q ) |0(g)) = 


moving-frame QRPA equations (local harmonic equations) 


(14) 


Sm ?)| [H M (q)iQ l {q)\ ~-^B ik {q)Pk(q) + \ 

1 k Z 

5 (</>(<?)I -Pm - Y c ij(q)Q j (q) -1 

X Zi 


dV 


_ k 


dq k 


m)) 


0,(15) 


f)V 
k 


,Y B ij(q)Q j (q ) 

j 


10 (g)) =0,(16) 


6 










where Hm (q) represents the Hamiltonian in the frame attached to the moving mean 
field, 


*»(!) =«- E 


(17) 


and is called ‘moving-frame Hamiltonian.’ The displacement operators Pi{q) and 
Cij(q) are defined by 


and 


A(g) 10(g)) =i -^i I0(«)) 


Cijiq) 


d 2 V ^ T kdV 

dq l dqi % i Qqk ’ 


(18) 


(19) 


respectively, with 



- b h ( ® Bh + i 

2 <9g* <90 


( 20 ) 


The double-commutator term in Eq. (11611 arises from the g-derivative of the infinites¬ 
imal generators Q l (q) and represents the curvatures of the collective submanifold. 

Solving these equations self-consistently, we can determine the microscopic ex¬ 
pressions of the infinitesimal generators Q l (q) and Pjfq) in bilinear forms of the 
quasiparticle creation and annihilation operators defined locally with respect to 
|0(g)). The collective Hamiltonian is given by 


with 


U(q,p) = (0(g,p)l-H'l0(g,p)) 

= v W) + \Yl 313 (vypiPj 

ij 


V(q)=H(q,p) 


B ij (q ) 


d 2 U 


dpidpj 


p =0 


( 21 ) 


( 22 ) 


where V(q) and (q) represent the collective potential and the reciprocals of col¬ 
lective inertial mass, respectively. They are functions of the collective coordinate 
q. Note that Eqs. (HU), (fT5|) . and (fl6l) reduce to the HFB and QRPA equations at 
equilibrium points, where dV/dq l = 0; namely, they are natural extensions of the 
HFB-QRPA equations to non-equilibrium states. 

Let us note the following points. 

• Difference from the constrained, HFB equations 

The moving-frame HFB equation (1141) resembles the constrained HFB equa¬ 
tion, but the infinitesimal generators Q l (q ) are here self-consistently deter¬ 
mined together with P t (q) as solutions of the moving-frame QRPA equations 
(fl5l) and (fT6l) at every point of the collective coordinate q. Thus, contrary to 
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constrained operators in the constrained HFB theory, their microscopic struc¬ 
tures change as functions of q. In other words, the optimal ‘constraining’ op¬ 
erators are locally determined at every point of q. The collective submanifold 
embedded in the TDHFB phase space is extracted in this way. The canonical 
quantization of the collective Hamiltonian described by a few collective vari¬ 
ables ( q,p ) is similar to the quantization of constrained system [30], but the 
‘constraints’ are here generated by the dynamics of the quantum many-body 
system under consideration. 

• Meaning of the term ‘adiabatic ’ 

It is used here in the meaning that we can solve the time-dependent variational 
equation (fT3l) in a good approximation by taking into account up to the second 
order in an expansion with respect to the collective momenta p. It is important 
to note that the effects of finite frequency of the LACM are taken into account 
through the moving-frame QRPA equations. No assumption is made like that 
the kinetic energy of LACM is much smaller than the lowest two-quasiparticle 
excitation energy at every point of q. 

• Physics of collective inertial mass 

The collective inertial mass represents the inertia of the many-body system 
against an infinitesimal change of the collective coordinate q during the time 
evolution of the mean field. It is a local quantity and varies as a function 
of q. As the single-particle-energy spectrum in the mean field changes as a 
function of q, level crossing at the Fermi energy successively occurs during 
the LACM. In the presence of the pairing correlation, the many-body system 
can easily rearrange the lowest-energy configurations at every value of q , i.e., 
the system can easily change q. The easiness/hardness of the configuration 
rearrangements at the level crossings determines the adiabaticity/diabaticity of 
the system. Since the inertia represents a property of the system trying to keep 
a definite configuration, we expect that the stronger the pairing correlation 
becomes, the smaller the collective inertial mass becomes m- It remains as 
an interesting subject to investigate how the self-consistent determination of 
the optimal directions of collective motion and the finite frequency ui(q) of the 
moving-frame QRPA modes affect the level-crossing dynamics of the superfluid 
nuclear systems. 

Consideration of gauge invariance 

In the QRPA, the ANG modes such as the number-fluctuation (pairing rota¬ 
tional) modes are decoupled from other normal modes and thereby the QRPA re¬ 
stores the gauge invariance (number conservation) broken in the HFB mean field 
[8]. It is desirable to keep such a merits of the QRPA beyond the small-amplitude 
approximation. Otherwise, spurious number-fluctuation modes would unexpectedly 
mix in the LACM of interest. We can take into account the gauge invariance in 
the following way [30]. Introducing the number-fluctuation variables ny't and the 
gauge angles conjugate to them, we write the TDHFB state vector (fl2l) in a 
more general form: 


\(f(q,p,<p,n)) = e 




\(/>(q,p,n)) 


(23) 


with 


(24) 


\<Kq,P,n)) = e^™ri\<j>(q)) 

and 

G(q,p,n) = Y^PiQ*(q) + Y ™ (t) © (t) (<?)• (25) 

i r=n,p 

Here N^) = fr(r) _ ]Sf( T ) are the number-fluctuation operators measured from 
N^ = {(f>{q)\ \(j)(q)), with the suffix r distinguishing protons and neutrons, 

and 0 ^ T )(g) infinitesimal generators for the pairing-rotation degrees of freedom. 
The state vector | (j>(q,p,n)} may be regarded as an intrinsic state to the pairing 
rotation. Using | (f)(q,p,ip,n)} in place of \4>(q,p)) in Eq. (fl3l) and expanding it in 
n as well as p up to the second order, we can determine &^(q) simultaneously 
with Q l (q ) and P l (q ) such that the moving-frame HFB+QRPA equations become 
invariant against the rotation of the gauge angle ipG) _ The gauge invariance of the 
resulting equations implies that we need to fix a gauge in numerical applications. 
A convenient procedure of the gauge fixing is discussed in Ref. m- 


3 Microscopic derivation of the 5D quadrupole 
collective Hamiltonian 

For collective submanifolds of two dimensions (2D) or more, large-scale numerical 
computation is needed to find fully self-consistent solutions of the ASCC equations. 
Then, a practical approximation scheme, called local QRPA (LQRPA) method, has 
been developed nn mm- This scheme may be regarded as a first step of iterative 
solution of Eqs. (fTTl) - (fl 6 l) . With use of it, we can easily derive the 5D collective 
Hamiltonian. We first derive the 2D collective Hamiltonian for vibrational motions 
corresponding to the ( 0 , 7 ) deformations, and then consider the three-dimensional 
(3D) rotational motions. 

First, we solve the constrained HFB equation 

<5 (0(0,7)1 # chfb (0, 7) 10(0,7)) = 0 (26) 


with 


^chfb(A 7) = ^-E A(t) (^^ (t) - E /bn(/?,7)^E 

r m=0,2 


where A (t) ( 0 , 7 ), // m (0,7) and D^ are the chemical potentials, the Lagrange mul¬ 
tipliers, and the quadrupole operators, respectively. The quadrupole deformation 
parameters ( 0 , 7 ) are defined by 


/3cos7 = = V ( 0 ( 0 , 7)1 10 (^, 7 )) , (27) 

-j=P sin 7 = r]D { 2 = 77(0(0, 7 )| D ( 22 ] |0(0,7)) , (28) 


where ?? is a scaling factor [TJ 01 ] . 


9 


Next, we solve 

8 (<K/3,7)I [#chfb(/ 3,7),Q*(A7)] - y-Pj(^,7) W.7)} = 0, (* = 1,2) (29) 

7)1 [£chfb(/?,7), tA(/?, 7)] - W,7)W,7) W,7)> =0- (i = 1,2) (30) 

These are the moving-frame QRPA equations without the curvature terms and called 
local QRPA (LQRPA) equations. Making a similarity transformation such that the 
collective masses corresponding to the collective coordinates ( 91 , 92 ) become unity, 
we can write the kinetic energy of vibrational motions as 

Fib = I E ( B ) 2 = i £ (<f) 2 (31) 

Z 1=1,2 ^ 1=1,2 

without loss of generality. Changes of the quadrupole deformation due to variations 
with respect to ( 91 , 92 ) are given by 

on(+) 

^ = E (rn = 0,2) (32) 

1 = 1,2 ^ 


Thus, the kinetic energy of vibrational motions is given in terms of time deriva¬ 
tives of the quadrupole deformation, 

T vib = ^Moo[^o } ] 2 + + \m 22 [D^] 2 , (33) 


where 


M mm i (p, 7 ) 


E 


dq l dq l 

ou 2m OU 2m' 


(34) 


It is straightforward to rewrite the above expression using the time derivatives of 
(/?, 7 ). Subsequently we solve the LQRPA equations for 3D rotational motions 
at every point on the (/3, 7 ) plane to obtain the inertial functions D}.{[3, 7 ) and 
the moments of inertia JkiPiZ) = 4/3 2 Z7 fc (/5, 7 ) sin 2 (7 — 2nk/3) determining the 
rotational energy ?rot- This step is the same as in Thouless and Valatin [45], except 
that the procedure is applied for non-equilibrium points of (/3, 7 ) as well as the 
equilibrium points in the potential energy surface V(/3, 7 ). 

After quantizing in curvilinear coordinates (so-called Pauli prescription) [46j, we 
obtain the quantized 5D quadrupole collective Hamiltonian, 


#coll — Tvib + Trot + V(/3, 7 ), 

whose vibrational kinetic-energy term takes the following form: 


TVib = 


+- 


- n 2 f 1 r & 


2 VWr [dp 


p c 



sin 3y 


d 


R 


R d_ 

IP 77 5/3 ^ 

d ' 


d 


-RZ, I V ^ sin3 T- D /37^ ) + 777 


57 


dp I $7 



(35) 


(36) 

(37) 
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where 


W = Dpp(/3, 7 )D 77 (fl, 7 )~D 2 7 (fl, 7 ), 
R = D 1 (fl, 7 )D 2 (fl, 7 )D 3 (fl, 7 ). 


(38) 

(39) 


If the inertial functions (Dpp,D 77 , D\, D 2 , D 3 ) are replaced by a common constant 
D and Dp 7 by 0, the above expression is reduced to 


Tv ib 


h 2 ( i 3 4 d 1 d . o d\ 

2D V/^d/T dp + P 2 sin3 7 Iff Sm ‘Hh() ' 


(40) 


Such a drastic approximation may be valid only for small-amplitude vibrations 
about a spherical HFB equilibrium. 

Writing the collective wave functions as 


v f'a7M(/3,7)^) 


(n\IMK) 


J2 QalldP^WMK), 

K=eve n 


21 + 1 


16vr 2 (l + 5 K o) 


\Pmk{&) + ( — YTm^k^I)] 


(41) 

(42) 


and solving the eigenvalue equation for vibrational wave functions 


T vib + V(P, 7 ) 


<S>aiK(P, 7) + E (IMK\ f rot | IMK') $ aIK '(P, 7) = E aI $ aIK (P, 7 ), (43) 

K'=eve n 


we obtain quantum spectra of quadrupole collective motion. The symmetries and 
boundary conditions of the vibrational wave functions are discussed in Ref. S3- 


4 Illustrative examples 

We here present some applications of the LQRPA method for deriving the 5D col¬ 
lective Hamiltonian. In the numerical examples below, the pairing-plus-quadrupole 
(P+Q) model Hamiltonian [38j (including the quadrupole-pairing interaction) is 
employed in solving the CHFB + LQRPA equations. The single-particle energies 
and the P+Q interaction strengths are determined such that the results of the 
Skyrme-HFB calculation for the ground states are best reproduced within the P+Q 
model (see Refs. [42, 49] for details). The LQRPA method is quite general and 
it can be used in conjunction with various Skyrme interactions or modern density 
functionals. A large-scale calculation is required, however, and such an application 
of the LQRPA method with realistic inter actions/functionals is a challenging future 
subject. A step toward this goal has recently been carried out for axially symmetric 
cases {50]. More examples can be found in [51J for 68 - ' 2 Se, [12] for 72 ’ 74 ’ 76 Kr, [II] 
for the 26 Mg region, |52j for 30 - 34 Mg, [50] for 58 - 68 Cr, [43] for 58 _ 66 Cr, and [53j for 

128-132^-g anc j 130—134g a 


Oblate-prolate shape coexistence and fluctuations in 72 Kr 
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Figure 1: Application of the LQRPA method to the oblate-prolate shape coexis¬ 

tence/fluctuation phenomenon in ,2 Kr (from Ref. |J2j)- (a) Collective potential V(/3, 7 ), 
(b) Ratios of the collective inertial masses Dpp(/3, 7 ) to the Inglis-Belyaev cranking 
masses. (c) Same as (b) but for D^(/3, 7 ). Vibrational wave functions squared, 
/3 4 \$ a iK{(3, y)| 2 , for (d) the 0 ] 1 " state, (e) the 2f state, (f) the 4f state, (g) the 0 ^ 
state, (h) the 2 \ state, and (i) the 4 2 state. For the /3 4 factor, see the text. 
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Figure 2: Excitation spectra and B(E 2) values calculated for 72 Kr by means of the 
CHB+LQRPA method (denoted by LQRPA) |42| and experimental data 'o4l [55] . For 
comparison, the results calculated using the Inglis-Belyaev cranking mass are also shown. 
Only B(E 2)fs larger than 1 Weisskopf unit are shown in units of e 2 fm 4 . 


The collective potential depicted in Fig. 1 exhibits two local minima. 

The oblate minimum is lower than the prolate minimum. This is expected from the 
deformed shell structure which gives rise to an oblate magic number at Z = N = 
36. This figure also shows that the valley runs in the triaxially deformed region 
and the barrier connecting the oblate and prolate minima is low. Accordingly, 
one may expect large-amplitude quantum shape fluctuations to occur along the 
triaxial valley. In fact, the vibrational wave function of the ground 0^ state peaks 
around the oblate potential minimum, but its tail extends to the prolate region. The 
vibrational wave function of the excited 0 ^ state consists of two components: one 
is a sharp peak on the oblate side and the other is a component spreading around 
the prolate region somewhat broadly. It is interesting to notice that, as the angular 
momentum increases, the localization of the vibrational wave functions in the (/?, 7 ) 
deformation plane develops; namely, the rotational effect plays an important role 
for the emergence of the shape-coexistence character. 

We note that not only the vibrational inertial masses shown in Fig. 1 but also 
the rotational inertial functions (D\, D 2 , and D 3 ) and the pairing gaps significantly 
change as functions of (/3, 7 ). Due to the time-odd contributions of the moving 
HFB self-consistent field, the collective inertial masses calculated with the LQRPA 
method are 20 — 50% larger than those evaluated with the Inglis-Belyaev cranking 
formula. Their ratios also change as functions of (fl, 7 ) [12]. As a consequence, as 
shown in Fig. 2, the excitation spectrum calculated with the LQRPA masses is in 
much better agreement with experimental data than that with the Inglis-Belyaev 
cranking masses. 

Quantum shape transitions and fluctuations in 30 > 32 > 34 Mg 
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This is a new region of quantum shape transition currently under live discussions 
toward understanding the nature of the quadrupole deformation in these neutron- 
rich isotopes as well as the mechanism of its emergence. 

Figure 3 shows the collective potentials V(f3, 7 ) and the vibrational wave func¬ 
tions squared, \^oIk{P, l )\ 2 1 of the 0^,2+, and states in 30 > 32 ’ 34 Mg. It is 
clearly seen that prolate deformation grows with the neutron number. The de¬ 
formed magic numbers, Z = 12 of protons and N = 20, 22, 24 at different values 
of (3 of neutrons [56] act cooperatively for the appearance of the prolate minima. 
Interestingly, the vibrational wave functions of the 2+ state are noticeably different 
from those of the 0j*~ state in 30 Mg and 32 Mg, while they are similar in 34 Mg. 

In Fig. 4, we display the probability density of finding a shape with a specific 
value of (3. 

P(i3) = f d7l^a7=o,K=o(^,7)| 2 |G(^,7)| 1/2 . (44) 

Note that the volume element with \G((3,'y)\ = 4/3 s kF(/3, -/)R(f3, 7 ) sin 2 37 is 
taken into account here. Let us first look at the ground 0^ states. The peak position 
moves toward a larger value of /3 in going from 30 Mg to 34 Mg. The distribution for 
32 Mg is much broader than those for 30 Mg and 34 Mg. Next, let us look at the 
excited 0^" states. In 30 Mg, the bump at /3 ~ 0.1 is much smaller than the major 
bump at /3 ~ 0.3. In this sense, we can regard the Oij - state of 30 Mg as a prolately 
deformed state. In the case of 32 Mg, the probability density exhibits a very broad 
distribution extending from the spherical to deformed regions up to (3 = 0.5 with a 
prominent peak at (3 ~ 0.4 and a node at /3 ~ 0.3. 

Thus, the shape coexistence picture that the deformed excited 0^“ state coexists 
with the spherical ground state approximately holds for 30 Mg. On the other hand, 
large-amplitude quadrupole-shape fluctuations dominate in both the ground and 
the excited 0 + states in 32 Mg, in contrast to the interpretation of gdeformed ground 
and spherical excited 0 + statesh [57] based on a naive picture of crossing between 
the spherical and deformed configurations. To test these theoretical predictions, an 
experimental search for the distorted rotational bands built on the excited 0^~ states 
in 30 Mg and 32 Mg is strongly desired. 


5 Some remarks on other approaches 

In Section 2, we reviewed the basics of a microscopic theory of LACM focusing on 
new developments achieved after 2000. In this section, we give short remarks on 
other approaches to LACM. Typical approaches developed by 1980 are described in 
detail in the textbook of Ring and Schuck [22], and achievements during 1980-2000 
are well summarized in the review by Do Dang, Klein, and Walet [55] . 

Constrained HFB + adiabatic perturbation 

Historically, the Inglis-Belyaev cranking masses derived from the adiabatic per¬ 
turbation theory [22] have been widely used in conjunction with phenomenological 
mean-field models, e.g. for the study of fission dynamics joTf . In recent years, it has 
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p p p 

Figure 3: Application of the LQRPA method to the low-lying states in 30 > 32 ’ 34 Mg (from 
Ref. [52]). (a) Collective potential V(j3, 7 ) for 30 Mg. The HFB equilibrium points are 
indicated by red circles, (b)-(d) Vibrational wave functions squared, \^oiIk(/3, 7 )| 2 , 
of the 0^,2+, and 0^ states in 30 Mg. Contour lines are drawn at every tenth part of 
the maximum value, (e)-(h) and (i)-(l): Same set of figures but for 32 Mg and 34 Mg, 
respectively. 
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Figure 4: (a) Probability densities integrated over 7 , P(/3) 

f d”/\$ aI=0 ,K=o(/3, 7 )| 2 |G'(/ 3 , 7 I 1 / 2 , for the ground Of states in 30 > 32 > 34 Mg plotted 
functions of f3. (b) Same as (a) but for the excited Of states. 


become possible to carry out such studies using self-consistent mean fields obtained 
by solving the constrained HFB equations [60]. The Inglis-Belyaev cranking masses 
have also been used for low-frequency quadrupole collective dynamics [ 6 U 1621[631163] , 
At present, a systematic investigation on low-lying quadrupole spectra is underway 
using the 5D collective Hamiltonian with Inglis-Belyaev cranking masses and the 
relativistic (covariant) density functionals [65, 1661 1B7, 68 i 69t 170117I J. 

A problem of the Inglis-Belyaev cranking formula is that time-odd mean-field 
effects are ignored and it underestimates the collective masses (inertial functions) 
[72). Moving mean fields possess time-odd components that change sign under time 
reversal operation, but the cranking formula ignores their effects on the collective 
masses. By taking into account such time-odd corrections to the cranking masses, 
one can better reproduce low-lying spectra [53 ]. 

Adiabatic TDHF theory 

Attempts to self-consistently derive the collective Hamiltonian using adiabatic 
approximation to time evolution of mean fields started in 1960 ’s mm • In these 
pioneer works, the collective quadrupole coordinates (/?, 7 ) were defined in terms of 
expectation values of the quadrupole operators and the 5D collective Hamiltonian 
was derived using the P+Q force model [08]. During 1970’s this approach was 
generalized to be applicable to any effective interaction. This advanced approach is 
called adiabatic TDHF (ATDHF) [71 ] [751 [76] . 

In the ATDHF theory, the density matrix p(t) is written in the following form 


as 
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(45) 

(46) 


and expanded as a power series with respect to \(i). 

p(t) = e ix{t) po(t)e~ ix ^ 

= Po(t)+i\x(t),po(t)] - ^L\(t), [x(t),Po(t)]\ + ... 

Baranger and Veneroni [74] suggested a possibility to introduce collective coordi¬ 
nates as parameters describing the time evolution of the density matrix po(t) and 
discussed an iterative procedures to solve the ATDHF equations. But this idea has 
not been realized until now. We note that the ATDHF does not reduce to the RPA 
in the small-amplitude limit if a few collective coordinates are introduced by hand. 
In fact it gives a collective mass different from the RPA m- 

The ATDHF theory developed by Villars [78] aims at self-consistent determi¬ 
nation of the optimum collective coordinates on the basis of the time-dependent 
variational principle. This approach, however, encountered a difficulty that we 
could not get unique solutions of its basic equations determining the collective path. 
This problem was later solved by treating the second-order terms of the momen¬ 
tum expansion in a self-consistent manner [371179]. It was shown that, when the 
number of collective coordinate is only one, a collective path maximally decoupled 
from non-collective degrees of freedom runs along a valley in the multi-dimensional 
potential-energy surface associated with the TDHF states. 

To describe low-frequency collective motions, it is necessary to take into account 
the pairing correlations. In other words, we need to develop the adiabatic TDHFB 
(ATDHFB) theory. This is one of the reasons why applications of the ATDHF 
have been restricted to collective phenomena where pairing correlations play minor 
roles such as low-energy collisions between spherical closed-shell nuclei [80 1 and giant 
resonances m- When large-amplitude shape fluctuations take place, single-particle 
level crossings often occur. To follow the adiabatic configuration across the level 
crossing points, the pairing correlation plays an essential role. Thus, an extension to 
ATDHFB is indispensable for the description of low-frequency collective excitations. 

In the past, Dobaczewski and Skalski [81] tried to develop the ATDHFB theory 
assuming the axially symmetric quadrupole deformation parameter /3 as the collec¬ 
tive coordinate. Quite recently, Li et al. [82] tried to derive the 5D quadrupole 
collective Hamiltonian on the basis of the ATDHFB. The extension of ATDHF to 
ATDHFB is not straightforward, however. This is because, as discussed in Section 
2, we need to decouple the number-fluctuation degrees of freedom from the LACM 
of interest and respect the gauge invariance with respect the pairing rotational an¬ 
gles. 

Boson expansion method 

Boson expansion method is well known as a useful microscopic method of describ¬ 
ing anharmonic (non-linear) vibrations going beyond the harmonic approximation 
of the QRPA. In this approach, we first construct a collective subspace spanned by 
many-phonon states of vibrational quanta (determined by the QRPA) in the huge¬ 
dimensional shell-model space, and then map these many-phonon states one-to-one 
to many-boson states in an ideal boson space. Anharmonic effects neglected in the 
QRPA are treated as higher-order terms in the power-series expansion with respect 
to the boson creation and annihilation operators. Starting from the QRPA about 
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a spherical shape, one can thus derive the 5D quadrupole collective Hamiltonian 
in a fully quantum mechanical manner. The boson expansion method has been 
successfully applied to low-energy quadrupole excitation spectra in a wide range of 
nuclei including those lying in regions of quantum phase transitions from spherical 
to deformed [82, M]. 

In the time-dependent mean-field picture, state vectors in the boson expansion 
method are written in terms of the creation and annihilation operators (r|,Tj) of 
the QRPA eigen-modes, or, equivalently, in terms of the collective coordinate and 
momentum operators ( Qi, Pi ), 


I <K*)> 


exp 


exp 


E (jh(t)Ti - r]*(t)T\ 
. i 

E 


qi 


l<M- 


(47) 

(48) 


With increasing amplitudes of the quadrupole shape vibration \r]i(t)\ (values of 
the collective coordinate qi(t)), anharmonic (non-linear) effects become stronger. 
Strong non-linear effects may eventually change even the microscopic structure of 
the collective operators ( Qi,Pi ) determined by the QRPA. In such situations, it is 
desirable to construct a theory that allows variations of microscopic structure of 
collective operators as functions of qi(t). It may be said that the SCC method has 
accomplished this task. 


Generator coordinate method 


The generator coordinate method (GCM) has been used for a wide variety of 
nuclear collective phenomena [85l 86, 18?] . Using the angular-momentum projector 
Pimk and the neutron(proton)-number projector Hy (Pz), we write the state vector 
as a superpositions of the projected mean-held states with different deformation 
parameters (j3, 7 ), 

| Vnzim) = [ mi E firziM i)PnPzPimk W, t)> • (49) 

Because the projection operators contain integrations, it has been a difficult task to 
carry out such high-dimensional numerical integrations in solving the Hill-Wheeler 
equation for the states |0(/3, 7 )) obtained by the constrained HFB method. In 
recent years, however, remarkable progress has been taking place, which makes it 
possible to carry out such large-scale numerical computations [88l \M I® ED M, 
193] , The HFB calculations using the density dependent effective interactions are 
better founded by density functional theory (DFT). Accordingly, the modern GCM 
calculation is referred to as ‘multi-reference DFT’ [88] . 

It is well known that one can derive a collective Schrodinger equation by making 
a gaussian overlap approximation (GOA) to the Hill-Wheeler equation [Ml [95l [96 ]. 
There is no guarantee, however, that dynamical effects associated with time-odd 
components of moving mean held are sufficiently taken into account in the collective 
masses (inertial functions) obtained through this procedure. In the case of center of 
mass motion, we need to use complex generator coordinates to obtain the correct 
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mass, implying that collective momenta conjugate to collective coordinates should 
also be treated as generator coordinates [221 [97]. 

A fundamental question is how to choose the optimal generator coordinates. 
With the variational principle, Holzwarth and Yukawa [98] proved that the mean- 
field states parametrized by a single optimal generator coordinate run along a valley 
of the collective potential energy surface. This line of investigation was further de¬ 
veloped [99] and greatly stimulated the challenge toward constructing a microscopic 
theory of LACM. In this connection, we note that conventional GCM calculations 
parametrized by a few real generator coordinates do not reduce to the (Q)RPA in the 
small-amplitude limit, differently from the case that all two-quasiparticle (particle- 
hole) degrees of freedom are treated as complex generator coordinates 1 00 . 

It is very important to distinguish the 5D collective Hamiltonian obtained by 
making use of the GOA to the GCM from that derived in Section 3 on the ba¬ 
sis of the ASCC method. In the latter, the canonical conjugate pairs of collective 
coordinate and momentum are self-consistently derived on the basis of the time- 
dependent variational principle. The canonical formulation enables us to adopt the 
standard canonical quantization procedure. Furthermore, effects of the time-odd 
components of the moving mean field are automatically taken into account in the 
collective masses (inertial functions). In view of the above points, it is highly desir¬ 
able to carry out a systematic comparison of collective inertial masses evaluated by 
different approximations including the ASCC, the ATDHFB, the GCM+GOA, and 
the adiabatic cranking methods for a better understanding of their physical impli¬ 
cations. In this connection, it is interesting to notice that the results of the recent 
GCM calculation for ' 6 Kr [92] . using the particle-number and angular-momentum 
projected basis (]49l) . are rather similar to those obtained by use of the 5D col¬ 
lective Hamiltonian with the Inglis-Belyaev cranking masses, except for an overall 
overestimation of the excitation energies by about 20%. 


6 Challenges for future 

As reviewed by Hyde and Wood [5], nature of low-lying excited 0 + states, system¬ 
atically found in recent experiments as well as those known from old days, is not 
well understood. It is thus quite challenging to apply, in a systematic ways, the 5D 
collective Hamiltonian approach to all of these data, from light to heavy and from 
stable to unstable nuclei, and explore the limit of its applicability. Recalling that 
the importance of the couplings between the quadrupole and pairing vibrations has 
been pointed out mmm, one of the basic questions is ‘under what situations 
we need to extend the 5D collective Hamiltonian to 7D by explicitly treating the 
proton and neutron pairing gaps as dynamical variables.’ 

Another interesting subject is to extend the collective Hamiltonian approach to 
a variety of collective phenomena, for example, those observed in rapidly rotating 
nuclei, heavy and super heavy nuclei, neutron-rich unstable nuclei, by taking into 
account the effects of rapid rotation and/or continuum, from the beginning in the 
single-particle (HFB) Hamiltonian. Macroscopic quantum tunnelings through self- 
consistently generated barriers, such as spontaneous fissions and deep sub-barrier 
fusions, are, needless to say, longstanding yet modern, fundamental subjects of 
nuclear structure physics. 
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It is a great challenge to develop the CHFB + LQRPA approach on the ba¬ 
sis of the time-dependent DFT. To efficiently solve the large-dimensional LQRPA 
equations containing density-dependent terms, the finite-amplitude method recently 
developed in Refs. [101. Il02i 11031 Il04j may be utilized. 
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